Methods for Predicting the Thermal Conductivity 
of Composite Systems: A Review 



R. C. PROGELHOF 
New Jersey Institute of Technology 
Newark, New Jersey 07102 

J. L. THRONE* 
Beloit Corporation 
Beloit, Wisconsin 

and 

R. R. RUETSCH 
New Jersey Institute of Technology 
Newark, New Jersey 07102 

The thermal conductivity of a solid or gas filled polymer is 
used in processing or end use application calculations. Numer- 
ous theoretical and empirical correlations are found in the 
literature. A carefulrevicw of these models indictes that no one 
correlation or technique accurately predicts the thermal con- 
ductivity of all types of composites. The investigation indi- 
cated that for solid filled, composites the Lewis and Nielsen 
equation fitted the experimental data best for the range of fillers 
tested. However, for a gas filled polymer, none of the theoreti- 
cal models proved adequate. The semi-empirical approach of 
Harding showed considerable merit. 



INTRODUCTION 

ry\ hroughout this review article a composite is defined 
-L as "a substance formed from one or more discrete 
components imbedded in a homogeneous, continuous, 
stationary medium". Specifically, our discussion is lim- 
ited to finite particles and/or gas filled voids imbedded in 
a medium. A group of discrete particles in intimate 
contact or particles in a stationary or moving air matrix 
are not considered to be composites. Early in the inves- 
tigation an attempt was made to classify all available 
references found in the literature into solid-solid and 
gas -solid compositions. However, it was apparent that 
the theoretical models were used by numerous authors 
for both types of systems. We, therefore, present herein 
many of the various theoretical and empirical models 
and provide a brief description of the basis for each in 
our literature review section. Discussion and comments 
are presented in the section entided General Discus- 
sion. 

LITERATURE REVIEW 

Listed alphabetically in the Bibliography are a group 
of papers associated with the energy transfer and/or gas 
diffusion in a filled composite. In this section are listed 
several theoretical or empirical models and a brief de- 
scription of their basis. Our literature review 
encompassed a considerable amount of material but in 
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no way should it be assumed complete. The following 
general nomenclature is used in the paper: 

k c = thermal conductivity of continuous phase 
k a — thermal conductivity of discrete phase 
k f = thermal conductivity of one of the discrete 
phases 

k e , k r> k v = thermal conductivity of composite 
0 = volume fraction of discrete phase 
x = parallel to filament axis 
y = perpendicular to filament axis 

Coefficients and other symbols unique to individual 
models will be identified when the model is analyzed. 
Series Model: 

k e = (l - 4>)K + (!) 

Parallel Model: 

L = 1 - j. + i. (2) 
k e k c kj 

Geometric Mean Model: 

k r = k d *k t "-*> (3) 

Behrens Theoretical Model (6) 

The theoretical results were based on a field solution 
to the equation of heat conduction for a composite mate- 
rial with orthorhombic symmetry. Solutions were pre- 
sented for elliptic filaments in a square lattice: 
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= i+(p-i)*-(p-i) : 



cr-f tl-^-(«r-D 
1 _ + 1)|--(P-1)^- 0 -) 



(4) 



= 1 + 



(p + X) -f- - (p - 1) W> - cr) 



(5) 



where V 



-4 



€ = 



Major axis of ellipsoid 
Focal distance 



For circular rods this reduces to: 

_ , f (P + D ± fezlM 1 
' " r L (p+D - (P~l)<h J 

And lor cubic symmetry of spheres: 

r te+2) + (p-iM 1 

*" - H (p+2) - (p-lW> J 
which reduces to Maxwell's equation. 



(6) 



(7) 



Bruggeman Theoretical Model (7) 

Using different assumptions than Maxwell for per- 
meability and field strength, Bruggeman derived the 
following implicit equation for dilute suspensions of 
spheres in a homogeneous medium: 



icit equation for dilute suspensions 
>mogeneous medium: 



(8) 



Budiansky Theoretical Model (9) 

Budiansky established an analogy between an elec- 
trostatic field and temperature distribution in a compos- 
ite consisting of a random mixture of TV isotropic con- 
stituents, at least (N - 1) of which are presumed to be 
roughly spherical in shape. The conductivity can be 
extracted from the following implicit equation: 



= i 



(9) 



3 k 



Cheng and Vachon Theoretical Model (12) 

Starting with Tsao's basic probablistic model, Cheng 
and Vachon assumed the discontinuous phase had a 
parabolic distribution. Based upon this assumed dis- 
tribution, Tsao's constants were evaluated and a closed 
form expressions for effective thermal conductivity of a 
two phase mixture as a function of discontinuous phase 
volume fraction were obtained: 



k e > k d 



1 



— tan" 1 



k e VC(k d -k e )[k e + B(k d -K)] 



b_ rci 

2 V]fc t .+ 



C(k a - K) 



K 

In 



B(k (1 - k c ) 
ka > k c 
1 



+ 



1 - B 



(10) 



VC(Jfc, - k d )[k c + B{k d - K)] 



B 



V(* c + B(k d - k c )) + ^VC{k r - * d ) 1 - B 



B 



+ 



(11) 



V*. + B(k d - k c ) - -rVC(k e - k d ) 



where for both equations: 

B = \fi$2 C = -4 y/2J(36) 

For three or more phase mixtures a discussion was pre- 
sented on several techniques that could be used to esti- 
mate effective thermal conductivities of the composite. 

Donea Variational Model (20) 

Using variational principles Donea obtained the up- 
per, k+, and lower, k~ t bounds lor the effective thermal 
conductivity of a composite material whichis statistically 
homogeneous and contains two phases each having iso- 
tropic and uniform thermal conductivity. 



k+ = <f>k~ + (1 - <f>)k v 
kk c 



k~ = 



k c + (1 - <t>)k 



(12) 
(13) 



Spheres: 



f U(P+2) + 2s(p-l)] 
(P+2) - (p-l)s 



(14) 



where p = k d lk c , s = {alb?, a = radius of inclusion, b = 
radius of the largest possible spherical shell surrounding 
the inclusion for random dispersion of spheres, andfc + = 

it" = k. 

Parallel Circular Fibers: 



Uto+1) +(p-1)jJ 
(p+1) - (p-l)s 



(15) 



where p =k d /k V) s = (a/b)' z , a = radius of fiber, and 6 = 
radius of largest hollow cylinder surrounding the fiber. 

These equations yield the lower bounds. For the case 
of no geometric constraints with p > 1 it has been shown 
that by reversing the roles of the continuous and discrete 
conductivities in the above equations form the upper 
bounds. 

Fricke Theoretical Model (24) 

Fricke extended Maxwell's analysis to homogeneous 
ellipsoidal particles that are far apart in a homogeneous 
medium: 
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' 1 + I F kd 



k = k 



1 + 0(F - 1) 



(16) 



where 



-I 



and 



The factor F represents the ratio of the overall tempera- 
ture gradients in the two phases and the//s the semi- 
principal axes of the ellipsoid. For a sphere f x = / 2 =fz 
and the equation reduces to Maxwells equation. 

Halpin-Tsai Theoretical Model (2) 

For filaments of uniform cross -sectional area arranged 
in parallel, the thermal conductivity parallel to the fila- 
ments was assumed to be (parallel model): 

k € = <f>k d + (1 - <t>)K (17) 

Using an analogy between in-plane field equations and 
boundary conditions to the transverse transport 
coefficient, the transverse thermal conductivity is: 



1 

' " 1 - T|0 



where 



hsL _ l 

H K 
v = r 1 

for plates of width a and thickness b: 

C-V5fc«(f) 

and for circular or square fibers: 

r = 1.0 

Hamilton and Crosser Semi-Theoretical Model (28, 29) 

Hamilton and Crosser define the thermal conductiv- 
ity of a two component mixture as: 



Mi - ♦) % 



, dr 



(18) 



The ratio of thermal gradients can be determined from 
the theoretical work of Maxwell and Fricke. 

dT 

nk 



dx 



L _ 



dT 

dx 



k d + (n- l)k t 



l _ t [ "*„ + (»- 0*, - (» - lfrfc - fcJ 1 ,, 9) 
6 ~ c L k„ + (n- l)h + *(* c - *„) J ^ 

where n is an empirical constant. For a spherical particle 
Maxwell showed analytically that n =3.0. For 
nonspherical particles Hamilton and Crosser correlated 
numerous sets of data with (0.58 < \fj < 1.0): 



n = 



where \(j is the sphericity. The sphericity was defined as 
the ratio of the surface area of a sphere, with a volume 
equal to that of the particle, to the surface area of the 
particles. 

In a similar manner the thermal conductivity of a 
multi-phase mixture was given as: 



k = k 



f, _ & »,(n, - 1) (K ~ k t )] 
I k ft, + (n, - lft -I 

r a Mi - M I 

L h k, + (n s - L)k,\ 



(20) 



Harding Semi-Empirical Model (32-36) 

Starting with an elementary conduction, convection, 
and radiation heat transfer analysis for low density, cel- 
lular foam, Harding obtained a semi-empirical relation- 
ship. 



K = 



+ kA + (k b - h)<M\ - F) (1 - 0.5IW)M 



c r v s 



U3 



T h + T C + 92 



1000 



(21) 



where E s = dimensionless efficiency fact for heat con- 
duction through "foam solids (empirical), k a = thermal 
conductivity of air, k b = thermal conductivity of blowing 
agent, R a = geometric ratio of cut surface exposed per 
unit foam volume, V y = volume of gas in average cell, M 
— mean mole fraction of blowing agent in a cell, and C r 
= coefficient defining the rate of radiant energy transfer 
through foam. 

Jefferson, Witzell and Sibbitt Theoretical Model (43) 

Jefferson, Witzell, and Sibbitt assuming the particles 
were spherical in nature devised a series-parallel equiv- 
alent system for a singular particle in a cube of continu- 
ous medium: 



= k e [l - 



77 



1 it [ (0.5 + n)k a k c ] 

4(1 + 2n) 2 J 4(1 + 2nf [ 0.5k c + nfcj 

(22) 



which substituted into Eq 18 results in: 
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Lewis and Nielsen Semi-Theoretical Model (53, 59-63) 
Lewis and Nielsen modified the Halpin-Tsai equation 
to include the effect of the shape of the particles and the 
orientation or type of packing for a two-phase system. 



w 



here 



** r \_l - B 4><Ij 1 



A = - 1 



(24) 



B = 



K 



<Jr = 1 + 



The values of A and <£ m for many geometric shapes and 
orientation are given in the following tables: 



Type of dispersed 
phase 



Direction of 
heat flow 



Cubes 
Spheres 

Aggregates of spheres 

Randomly oriented rods 

Aspect ratio = 2 
Randomly oriented rods 

Aspect ratio = 4 
Randomly oriented rods 

Aspect ratio = 6 
Randomly oriented rods 

Aspect ratio = 10 
Randomly oriented rods 

Aspect ratio = 15 
Uniaxially oriented fibers 
Uniaxially oriented fibers 



Any 
Any 
Any 

Any 

Any 

Any 

Any 

Any 



2.0 

1.50 

2.50 

+. 
1.58 

2.08 

2.8 

4.93 

8.38 



-1 



Parallel to fibers 2L/D 
Perpendicular to 
fibers 0.5 



Shape of Particle Type of Packing 



4> m 



Spheres 
Spheres 
Spheres 
Spheres 
Spheres 
Spheres 
Rods or fibers 
Rods or fibers 
Rods or fibers 
Rods or fibers 



Hexagonal close 0.7405 

Face centered cubic 0.7405 

Body centered cubic 0.60 

Simple cubic 0.524 

Random close 0.637 

Random close 0.601 

Uniaxial hexagonal close 0.907 

Uniaxial simple cubic 0.785 

Uniaxial random 0.82 

Three dimensional random 0.52 



The values of A listed in the table are for elastic moduli. 
In general, the theoretical values for conductivity are 
slightly greater. 



Loeb Theoretical Model (54) 

Loeb derived a relationship for the effective thermal 
conductivity of a porous solid as a function of the cross- 
sectional and longitudinal pore fractions and surface 
emissivity. 

*• ° - ; w?v - 5 )1 (25) 

where <f> = cross-sectional porosity, fa = longitudinal 
porosity, y = geometrical factor, e ==_pore emissivity, 
cr = St^fan-Boltzman constant, and AT = temperature 
difference. 

The Loeb model attempts to compensate for radiation 
in the direction parallel to heat flow. For an orderly 
array of pores the above model was incorporated with a 
series-parallel analysis. 

Maxwell Theoretical Model (56) 

Maxwell using potential theory obtained an "exact" 
solution for the conductivity of randomly distributed 
and non-interacting homogeneous spheres in a 
homogeneous continuous medium: 

_ r fc, + 2h + Ufa - y 1 {26) 

Norton Empirical Model (64) 

Norton's techniques were based upon the following 
assumption. The polymer matrix which makes up the 
foam network is so complex that it would be impossible 
to model the conductive and radiative energy transfer 
within the foam. The thermal conductivity of a polymer 
foam with a known gas content is experimentally mea- 
sured. Subtracting the thermal conductivity of the gas 
assuming it had the same total volume would give the 
"effective" conductivity of the polymer for that particu- 
lar structure and density. 



^effvolvmer ~ ^measured 



- k 



gas 



The thermal conductivity of the foam with a different gas 
is then assumed to be 



= k 



+ k, 



■foam - KefI tJO i vme r T *U a * 



(27) 



Peterson and Hermans Semi-Theoretical Model (67) 

Peterson and Hermans derived a general theory for 
the dielectric constants of a dilute suspension of spheres 
in a continuum. By direct analogy between heat transfer 
and electrostatics. 



K = KW + *y<t> + V 



(28) 



where 



k d — k c 

7 ~~ 2k c + Jfc d 
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Progelhof and Throne Empirical Model (70) 

By curve fitting the experimental data for several 
foamed plastics they found that the following empirical 
equation fit the data best: 

k e = k d [l +A (fJ] (29) 

where 

and B is an empirical constant. 

Ratcliffe Empirical Model (71) 
Starting with a geometric mean correlation: 

Ratcliffe analyzed approximately 150 sets of data and 
determined the variation of the individual data points 
with the correlation. Based on a mean variation effect he 
prepared a chart defining probable limits (upper and 
lower) to the effective thermal conductivity of the two 
phase medium. Based upon these results he presented a 
nomograph for predicting effective thermal conductiv- 
ity. 

Russell Theoretical Model (75) 

Assuming the pores are cubes of the same size and the 
isothermal lines are planes, Russell obtained the con- 
ductivity using a series parallel network: 



k = k 



4* + f - (i - 4> m ) 



- <f> + f- (1 



+ 0-^ 3 )j 



(30) 



Springer and Tsai Semi-Theoretical Model (78) 

The thermal conductivity of a composite with uni- 
directional filaments in the direction along the filaments 
(parallel model) 

In the direction normal to the filament an analogy be- 
tween the response of an unidirectional composite to 
longitudinal shear loading and heat transfer results in: 



• (2a - h) + 



] 



where s = maximum transverse dimension of filament, 
h = width of the filament, and a, b = dimensions of 
elemental volume in which filament is mounted. 
For a square filament in a square array: 



(31) 



and cylindrical filaments in a square array: 
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Tan 



VE5 



i + 




(32) 



w 



here 



Topper Theoretical Model (81) 

Using a series-parallel energy model for a system of 
cubic voids in orderly array, the equivalent thermal 
conductivity was found to be: 

eft*/ 3 



K'i k c + 



-#4 



k d <t>*« + K(i 

For spheres distributed in a simple cubic lattice: 



(33) 



- m° 



Tan 



1/3 



r fc c _ 


'3$ 


2/3" 


1/2 




An 







(34) 



Radiant energy transfer has been neglected in this 
model. 

Tsao Theoretical Model (83) 

Tsao's model is based on determining a one dimen- 
sional porosity, P u and the mean, //,, and standard devia- 
tion, O", which relate P x to the "two dimensional" or 
planar porosity, P 2 . The bulk porosity, <f>, is the product 
of F x X F 2 . Using a probablistic model Tsao developed a 
relationship between P x and P 2 . Using a parallel conduc- 
tive model the equivalent thermal conductivity is given 
by the following integral equation: 



JL = 



f 1 dF* 



(35) 



To evaluate this integral the relationship between >x, 
and cr must be known. 

GENERAL DISCUSSION 

The general discussion on the appliability of the 
theoretical/empirical models presented in the preced- 
ing section has been divided into two parts. The first part 
deals with solid-solid compositions and the second with 
gas-solid composites. To evaluate the various models 
presented in the previous section, it is necessary to 
obtain reliable experimental data. In order that data be 
acceptable, they must satisfy several important criteria. 
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First, the experimenter must have described in detail a 
suitable technique or method for preparing test speci- 
mens. Second, the data must be of sufficient size and 
range of the concentrations to be useful for evaluation at 
both ends of the discontinuous phase concentration and 
must have included values of the continuous phase 
alone. And most important, a reliable, reproducible 
experimental technique and test set-up mast have been 
used to measure the thermal conductivity of the mate- 
rials. Much of the available experimental data given in 
the literature unfortunately do not meet any or all of 
these criteria. This does not mean that the experimental 
values not used in this analysis are incorrect, only that 
they are not acceptable for the purpose of evaluating the 
accuracy of the various theoretical models. 

Solid-Solid Compositions 

For a two component mixture the authors have 
selected the data of Sundstrom & Lee (80) for 62-88 
micron glass spheres and 62-125 micron magnesium 
oxide powder in polyethylene. These experimental 
data, shown in Figs, 1 and 2, were chosen by the authors 
since the range of filler to matrix conductivity ratio 
varied from 3 for glass to 164 for magnesium oxide. 
Measurements were made for filler contents of from 
zero to 30 percent in 5 percent increments. Also shown 
in the figures are the theoretically predicted variations 
of thermal conductivity by the following models: 

1. Bruggemen: Eq 8 

2. Cheng & Vachon: Eqs 10 and 11 

3. Jefferson, Witzell 6r Sibbitt: Eqs 22 and 23 

4. Lewis & Nielsen: Eq 24 

5. Maxwell: Eq 26 

6. Russell: Eq 30 

The theoretical values predicted by the Lewis and 
Nielsen model are based on <p m = 0.637, A oUm — 1.5, 
and A M = 3.0. Many of the other models reviewed in 
the Literature Review section were not plotted since 
these models reduce to the Maxwellian model when die 
particles are spherical in shape. 

After careful examination of curves for all fillers (79, 
80) over the complete range of fillers and additional data 
of a limited scope (11, 28, 71, 76), it is the authors 
conclusion that the Lewis and Nielson correlation fits 
the data best for the range of fillers tested. 

For a two component composite of highly non- 
spherical particles, such as glass filaments, we have been 
unable to obtain reliable data for direct comparison with 
the various theoretical models. The majority of models 
found in the literature for circular or elliptical fibers 
were derived for specific geometries and spacings of the 
filaments. The correlation of Lewis and Nielsen does 
account for filament size, orientation, etc. However, the 
data we evaluated as reported by Ashton et at. (2) and 
Segal (76) did not accurately describe the specimen 
preparation for die thermal conductivity measurement 
system. Nielsen presented an analysis of the Ashton et 
at. data and showed that for the possible range of vari- 
ables of the physical system his correlation would repre- 
sent the data. 



12 




Volume percent glass 



Fig. i. Glass filled polyethylene. A. Bruggman, B. Cheng 6 
Vachon, C.Jefferson et al., D. Maxwell, E. Lewis 6 Nielsen, F. 
Peterson & Hermans, and G. Russell. 




10 20 30 

Volume percent MgO 



Fig. 2. Magnesium oxide filled polyethylene. A. Bruggman, B. 
Cheng <b Vachon, C.Jefferson et aL, D. Maxwell, E. Lewis 6 
Nielsen, t\ Peterson 6- Hermans, and G. Russell. 

A more difficult and perplexing problem is the estima- 
tion of the thermal conductivity of a composite of a 
mixture of fillers. This system is encountered in many 
practical applications. Most likely, one of the fillers 
would be rod shaped in geometry for strengthening. 
Two correlations reviewed by the authors for multi- 
phase composites are those of Budiansky (9) and Hamil- 
ton (29). Budiansky \s model assumes that the particles 
are near spherical in shape. Hamilton's model uses a 
sphericity function to account for non-spherical parti- 
cles. However, there is no orientation effect. Limiting 
discussion to near spherical particles, specifically to- 
glass spheres and MgO powder, the authors predict the 
thermal conductivity of filled polyethylene by Budian- 
sky's, Hamilton's and Lewis & Nielsen's models. In the 
Hamilton model the sphericity of the glass spheres is 
assumed to be 1.0, and the MgO powder is assumed to 
be 0.625. Since the Lewis and Nielsen correlation is for a 
two component composite mixture, the three compo- 
nent mixture was estimated by first determining the 
conductivity of the polymer and one filler and then the 
conductivity of the mixture with the third filler. Results 
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were obtained for the fillers being added in different 
sequences. The results are listed in Table 1 . Based upon 
these results, it is apparent that the order of filler addi- 
tion when using the Lewis and Nielsen model has little 
effect. The results obtained by Budiansky's model are in 
fair agreement with the Lewis and Nielsen model but 
both significantly underestimate the composite conduc- 
tivity. Owing to the lack of accurate data, the authors are 
reluctant to recommend any correlation. 

In predicting the thermal conductivity of any two or 
multi phase mixture by one of the theoretical or empiri- 
cal relations discussed, hi this paper, it is assumed that 
good adhesion between the polymer and the filler is 
present. If there is a thin gas layer between the filler and 
the polymer, the "effective" thermal conductivity of the 
particle could be considerably smaller than the thermal 
conductivity of the pure filler phase. Micro photographs 
(17, 73) of fractured glass fiber reinforced material indi- 
cates that with poor wetting this factor would be apprec- 

* 3 ^ e * Gas- Solid Compositions 

In the preceding section the effect of solid additives 
on the thermal conductivity of a continuous matrix was 
discussed. The only energy transfer process that occurs 
within the composite, if complete wetting between the 
filler and the continuous media is assumed, is conduc- 
tion. However, the energy transfer process in a gas filled 
composite is much more complex. Shown in Fig. 3 is an 
idealized model of a gas filled composite filling the space 
between two parallel planes which are at temperatures 
T, and T 2 . respectively. The insert is of a single bubble, 
assumed to be spherical in shape. First assume that the 
continuous phase is opaque to electromagnetic radia- 
tion. This important assumption will be discussed 
further in our analysis. The only mode of energy transfer 
in the continuous phase is conduction. If the gas bubbles 
are relatively small and spaced far apart, the isotherms 
or lines of constant temperature for steady state condi- 
tion will be parallel to planes 1 and 2. These isotherms 
are shown as dotted lines. Thus, the temperature at any 
point on the surface of the bubbles or in the polymer 
assuming the thermal conductivity is constant is given 
by: 

Table 1. Thermal Conductivity 

Glass content, 
percent 



MgO content, 
percent 


0 




10 




20 


0 


8.01 


1. 


8.98 


1. 


10.13 




1. 11.59 


1a. 


13.1 


1a. 


15.0 


10% 




b. 


13.0 


b. 


14.6 




2. 


9.4 


2. 


10.0 






3. 


13.1 


3. 


15.1 




1. 16.50 


1a. 


19.0 


1a. 


22.1 


20% 




b. 


18.4 


b. 


20.6 




2. 


10.28 


2. 


10.92 






3. 


22.6 


3. 


26.5 



© 




Isotherms 



■H 



Fig. 3. Gas filled polymer. 



r = r t + T \ Tx * 

where x is the distance from plane 1. With a linear 
temperature gradient across the bubble, energy will be 
transferred by conduction and/or convection inside the 
bubble and by radiation on the inside surface of the 
bubble. If the gas molecules within the spherical bubble 
transfer energy only by molecular action, the rate of 
energy transfer is governed by Fourier's law: 

dT 



A ~ 



dx 



1. Lewis & Nielsen 

a. (Polymer & Glass) 4 MgO 

b. (Polymer & MgO + Glass) 

2. Hamilton & Crosier 

3. Budlantky 



If, however, the energy level (and thus the temperature 
difference) across the bubble is large, energy transfer 
will occur by molecular motion as well as by molecular 
transport (8, 14). The rate of energy transfer is governed 
by Newton's law of cooling: 

j^=h (T„ - T h ) 

where h is the convective conductance and T n and T b are 
the temperature of the surface and the fluid respec- 
tively. The molecular motion will result in a microscopic- 
velocity field developing inside of the bubble (13), Thus, 
the rate of energy transfer will be significantly greater 
than that which would occur with pure molecular con- 
duction. The development of the internal velocity field 
which enhances the rate of heat transfer is dependent 
upon many factors. These include the bubble diameter, 
temperature gradient, orientation, gas properties, etc. 
Techniques to correlate the effect of each of these pa- 
rameters on the rate of energy transfer usually rely on 
non-dimensional numbers. The common correlation for 
free convection effects is (7): 
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Nusselt Number Grashof-Prandtl Number 

where: p = gas density, £ = coefficient of thermal 
expansion of the gas, fj, = absolute viscosity of the gas, d 
= bubble diameter, and AT„ = temperature difference 
across the bubble. 

Thus, for a given fluid the significant parameters are 
d 3 and AT b . However, the temperature drop across the 
bubble is a function of the overall temperature drop' 
across the foam sheet: 



Therefore, 



¥-/(• 



L 

It is apparent that the convective effects are a function of 
the temperature, gradient, direction of energy How 
(horizontal or vertical) and the bubble diameter to the 
fourth power. The only criterion commonly found in the 
literature is the bubble diameter. Convective effects are 
largest for large diameter bubbles. Effective cell or bub- 
ble size for relatively low density foams can be found in 
Table 2 of Guenthers (26) paper. From this table it is 
evident that cell size is inversely proportionally to foam 
density. Thus, the convective effects will be greatest in 
low density foams. 

An estimate of the relative magnitude of the energy 
transfer within a bubble by radiation as compared to that 
by pure conduction can be made by considering the 
energy transfer between two parallel planes spaced a 
distance apart. This geometry is used instead of a sphere 
in that it can be analyzed theoretically. The planes are 
assumed to be perfect absorbers and at temperatures T a 
and T b . If T a and T b only differ by a temperature differ- 
ence e, T a can be rewritten as: 

The rate of energy transfer between the two surfaces is 
described by: 

conduction: -~ = ^- (T a — T b ) 

radiation: -2- = or(T* - T b ) 
A 

where o* is the Stephan-Boltzman constant. Since the 
temperature difference across the surfaces is e the con- 
ductances for the two modes of heat transfer are: 

g -i* 

^conduction J 



Gcowcion = 0.6872 Xl0- 8 n€ 

Shown in Fig. 4 are comparisons using air as the 
encapsulated gas for various ranges of bubble sizes in 
many foam applications. Even though these results were 
derived for a thin air gap between two planes, the rela- 
tive effects would be similar for a spherical inclusion. 
Assuming no internal convective currents, the results 



clearly indicate that radiation can be a significant portion 
of the total energy transfer within a cell when the cell 
size is large, d - 0.040 in. , and the temperature is high, 
T > 200°F. As was shown in the section on convective 
effects, the greatest effect will occur in low density 
foams. 

If the energy transfer process in a cell is treated as the 
sum of convection and radiation, there will be a specific 
cell size at which the effective thermal conductivity of 
the foam will reach a minimum and then increase as the 
foam density is further decreased. This effect is clearly 
seen in the data of Baxter and Jones (5) for polystyrene, 
which has been replotted in Fig. 5. The minimum ther- 
mal conductivity, 0.0173 Btu/hr ft°F, occurs at foam 
densities of from 2-8 Ib/cu ft. It is also apparent that 
foams will exhibit significantly different thermal proper- 
ties depending upon the direction of foam rise. In effect 
the cells will not be true spheres but will be ellipsoidal in 
shape. The convective effects in the direction of the 
major axis of the inclusion will be greater than that in the 
direction of the minor axis. Since the direction of rise is 
the same as the major axis the effective thermal conduc- 
tivity in the direction of rise will be greater than that of 
the transverse direction. In addition, if the foam is an 
open cell geometry, the air may migrate from cell to cell 
forming a natural convection pattern within the foam 
body. 
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Fig. 4. Conduction/radiation effects. 
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Fig. 5. Thermal conductivity of polystyrene foam. 
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One of the original assumptions made in the analysis 
of the energy transfer process in a gas filled composite 
was that the continuous phase does not transmit elec- 
tromagnetic radiation. However, this assumption is not 
valid for polymers over all conditions. The characteristic 
property that identifies the degree of transparency of a 
material is the absorption coefficient, a, which is de- 
fined in the following manner: 

dZ(X) = - a(k) I(k) dx 

where: I = incident energy beam, A = wavelength, and 
x = distance. As a decreases in value, the transmittance 
of radiant energy through a sheet of given thickness 
increases- Tf a is zero at a given wave length, the sheet is 
totally transparent, regardless of sheet thickness. Nor- 
mally for a given material the absorption coefficient is a 
highly irregular function of wavelength. This is seen in 
Fig. 6. Also shown in the figure is an approximation of 
the irregular curve by a step function curve. The data 
presented in Fig. 6 were obtained from sheet speci- 
mens. Results for thin films can be significantly different 
than those for a sheet. The apparent discrepancy occurs 
when the wave length of the incident radiation is of the 
same order of magnitude as the film thickness. 
Additional unpublished data for nylon-6 and 
polyethylene sheets indicate that crystalline polymers 
are stronger absorbers than noncrystalline polymers. 
However, most polymeric materials do transmit some 
radiant energy. It is necessary to determine the rate of 
transmission for various materials, however. 

The simplest model to describe the transmission of 
radiant energy through an absorbing-scattering media is 
the "two flux" or one-dimensional model. A comprehen- 
sive discussion of this technique can be found in refer- 
ences (30, 69, 74). Estimating the rate of energy transfer 
in a scattering absorbing medium is extremely complex 
and is a function of the type of inclusion, polymer, 
radiation spectrum, etc. It can be that there is no scatter- 
ing and if the temperatures of the polymer and the 
source are low, (e.g., less than 200°F) and if the foam cell 
walls have absorption characteristics shewn in Fig. 6. 
The results of Harding and James (31) indicate that a low 
density polystyrene foam may have a wall thickness of 
only 40 molecules thick. This thickness apparently is a 
lower limit to differentiating between a film and a sheet. 
This conclusion is supported by the data of Chute, 
Walker and DeWitz on the radiant heating of low den- 
sity polystyrene foam indicating that their material was 
relatively opaque to long wavelength radiation but con- 
siderably transparent to short or high temperature radi- 
ation. 

Most of the theoretical models listed in the beginning 
of the report are based on a specific thermal conductivity 
of the discrete phase of the composite. From our discus- 
sion it is apparent that this idealization is not representa- 
tive of many low density foams and application of these 
models would result in erroneous predictions. Each 
polymer, gas, application, etc. must be considered be- 
fore attempting to estimate the thermal conductivity of 
the foam. For this reason the approach of Harding (32- 
36), based on a semi-empirical approach, shows consid- 
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Fig. 6. Absorption characteristics of polystyrene. 

erable merit. However, for each polymer, configura- 
tion, and environmental condition the empirical con- 
stants should be evaluated. For example, if empirical 
coefficients were measured for an effective thermal con- 
ductivity of polystyrene in a room temperature 
environment, the application of these coefficients to a 
radiant heating system would also result in erroneous 
energy transfer rates. 

In our discussions we have based our assumptions that 
the properties of the encapsulated gas are known. How- 
ever, when foams are blown with gases other than air 
there is an aging effect (16, 18, 19, .36, 42, 52). This is 
due to the diffusion of blowing gas out of the foam and air 
into the foam. Recently models have been proposed to 
account for the diffusion process, The models usually 
assume that the cells are cubic and the gas diffuses 
through a series of sheets. This model appears to be 
satisfactory for low density foams but may not be applic- 
able for higher density regions, such as the skin of a 
structural foam. 
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